Random Matrix Theory of a Chaotic Andreev Quantum Dot 
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A new universality class distinct from the standard Wigner-Dyson ones is identified. This class 
is realized by putting a metallic quantum dot in contact with a superconductor, while applying 
a magnetic field so as to make the pairing field effectively vanish on average. A random-matrix 
description of the spectral and transport properties of such a quantum dot is proposed. The weak- 
localization correction to the tunnel conductance is nonzero and results from the depletion of the 
density of states due to the coupling with the superconductor. Semiclassically, the depletion is 
caused by a singular mode of phase-coherent long-range propagation of particles and holes. 
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Wigner-Dyson level statistics is found in physical sys- 
tems as diverse as highly excited molecules, atoms and 
nuclei, mesoscopic systems in the ballistic or diffusive 
regime, and chaotic Hamiltonian systems such as the sta- 
dium or Sinai's billiard. The reason for the ubiquity and 
universality of Wigner-Dyson statistics is the relation of 
the Gaussian ensembles |IJ to attractive fixed points of 
the renormalization group flow for an effective field the- 
ory (nonlinear a model) || . Depending on whether time 
reversal and/or spin rotation invariance is broken or not, 
the relevant ensemble has orthogonal, unitary, or sym- 
plcctic symmetry. 

Although the Wigner-Dyson ensembles are the generic 
ones, new universality classes may arise when additional 
symmetries or constraints are imposed. One example of 
this is provided by systems with A-B sublattice structure 
or chiral symmetry ||. Another example, presented in 
this letter, are metallic systems in contact with a super- 
conductor. Our considerations were inspired in part by 
work [Q on Anderson localization of normal excitations 
in a dirty superconductor or a superconducting glass. 




FIG. 1. Metallic quantum dot (N) in contact with a super- 
conductor (S). A tunnel barrier (T) separates the dot from 
a normal-metal lead (L) attached for the purpose of making 
current and voltage measurements. 

The system we have in mind is depicted in Fig. 1: a 
metallic (normal-conducting) quantum dot (N) with the 
shape of, say, a stadium billiard, is surrounded by a su- 



perconductor (S) and contacted by a normal-metal lead 
(L) via a tunnel barrier (T). A Schottky (or potential) 
barrier forms at the NS-interface. The quantum dot is 
pierced by a weak magnetic flux of the order of one or 
several flux quanta. The temperature is so low that the 
phase-coherence length exceeds the system size by far. It 
will help our argument if the quantum dot is not perfectly 
shaped but contains some defects and/or impurities. 
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FIG. 2. (Almost) periodic orbit involving several normal 
reflections (N) and one Andreev reflection (A). An electron 
runs through the closed loop twice. Before reaching A it is a 
particle and afterwards a hole. 

The unique feature that distinguishes our quantum dot 
from conventional mesoscopic systems is the process of 
Andreev reflection [pi, by which an electron incident on 
the NS-interface is retroreflected as a hole, and vice versa. 
The drastic consequences of this process for the thermo- 
dynamics and the transport of the quantum dot can be 
anticipated from a semiclassical argument. Consider the 
periodic orbit of Fig. 2: an electron starts out from point 
A, undergoes several normal reflections off the Schottky 
barrier and eventually returns to A, with the final veloc- 
ity being roughly the negative of the initial one. At the 
point of return, the electron gets converted into a hole 
by Andreev reflection. Because the magnetic field is too 
weak to cause any significant bending of classical trajec- 
tories, the hole then simply tracks down the path laid 
out by the electron. The charge of the hole is opposite to 
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that of the electron, so the magnetic phase accumulated 
along the total path is zero. Moreover, since a hole is 
not just the charge conjugate but also the time reverse of 
an electron, the dynamical phases cancel, too, provided 
that electron and hole are at the same energy (the Fermi 
energy). The two Andreev reflections add up to a total 
phase shift of n. Thus, the periodic orbit of Fig. 2 con- 
tributes to the periodic-orbit sum for the density of states 
with the negative sign. Orbits of this type reduce the 
mean density of states ("Weyl term") and are expected 
to put our "Andreev quantum dot" in a universality class 
distinct from the standard Wigner-Dyson ones. We shall 
identify this universality class in the sequel. 

In a microscopic mean- field treatment, we would start 
from the Bogoliubov-deGennes (BdG) Hamiltonian 

Tt=f H ° A \ 

where Hq = (p — eA) 2 /2m + V(x) — /i is a Hamiltonian 
for "particles", — Hq is the corresponding Hamiltonian 
for "holes" , and the pairing field A(x), whose magnitude 
rises from zero inside the quantum dot to a nonzero value 
in the superconductor, converts particles into holes. The 
potential U(x) includes the Schottky barrier. Energy is 
measured relative to the chemical potential //. 

Universality classes are characterized by their symme- 
tries. Notice therefore that 7i obeys the relation 

n = -cn T c-\ c=(\lY (i) 

This symmetry originates from the electron's being a 
spin- 1/2 fermion and from spin-rotation invariance 
The transformation H h-> CJ^C^ 1 will be called CT- 
conjugation as it combines time reversal (Ti i— > Ti* = 7Y T ) 
with a kind of charge conjugation (Ti i— > CTiC^ 1 ). 

By design, the classical motion of particles and holes 
in the billiard-shaped dot is chaotic and fills the avail- 
able phase space ergodically. Now observe that every 
time a particle or hole is Andreev-reflected from the NS- 
interface, its wavefunction acquires an extra phase deter- 
mined by the superconducting order parameter. In this 
context it is important that the applied magnetic field 
is screened by a supercurrent circulating along the NS- 
interface inside the superconductor. The supercurrent 
flow, in turn, is concomitant with a spatial variation of 
the phase 4> of the order parameter |Q. As a result of 
this and the chaotic dynamics, the extra phase picked up 
during Andreev reflection varies randomly along a typical 
semiclassical trajectory. So everything is quite random, 
and we expect some kind of random-matrix theory to ap- 
ply. The question is now: what random-matrix theory? 

Because the presence of the magnetic field makes the 
pairing field experienced by particles and holes vanish 
on average, A can be modelled by a stochastic vari- 
able with zero mean. Moreover, since the system has 



been designed to be chaotic, there exist no integrals of 
motion except for energy. The only symmetry (apart 
from hermiticity of the Hamiltonian) of relevance for 
the long-time ]8j or ergodic limit we shall consider, is 
the CT-oddness (Q). Experience with similar problems 
then tells us that we can model the ergodic limit by 
a Gaussian, or maximum-entropy, ensemble with prob- 
ability density exp(— TrH 2 /2v 2 )dH subject to the con- 
straint (|l|). This implies that, in any orthonormal basis 
of states ip a (a — 1,1, N, N) with CT-conjugate ba- 
sis ipa = Cip*, the variances of the random Hamiltonian 
matrix elements are given by the correlation law 

(HabHcd) = v 2 {S ad S bc - 5 a c6 b3 ). (2) 

To complete the definition of our random-matrix model, 
we add to Ti a term — iT which accounts for the coupling 
to the normal-metal lead and will be specified later. 

As a first step, let us close off the contact with the lead 
(r = 0). What can we say about the spectral statistics, 
the central characteristic of the ergodic isolated quantum 
dot? If ipk is an eigenstate of Ti with eigenvalue +Ek 
then, by (0), so is Cip% with eigenvalue — E}.. Thus, there 
is an exact pairing between positive and negative eigen- 
values. By diagonalizing Ti and computing the Jacobian 
of the transformation to diagonal form, we obtain the 
(unnormalized) joint probability density of the positive 
eigenvalues Ek (k = 1, N); 

N 

P(E U E N ) = Y[(E? - E 2 ) 2 J] E\ e- E ^, (3) 

i<j fe=l 

which is manifestly invariant under Ek — > —Ek 0. 

To calculate the spectral statistics, it is convenient 
to view (H) as a Gaussian Unitary Ensemble (GUE) of 
2N levels Ei,Ei,...,E^,Ej^ with the mirror constraint 
Ej. = —Ek- The correlation functions of the GUE are 
known to coincide with those of a one-dimensional gas 
of free fermions in the large- N limit . Now, when a 
fermion (i.e. an energy level) gets close to E = 0, so 
does its mirror image. Because Fermi statistics makes 
the wavefunction vanish as two fermions approach, the 
constraint Ej. = —Ek amounts to hard wall boundary 
conditions at E = 0. Hence, we can compute the eigen- 
value density and its correlations for (|J) as the particle 
density and its correlations for a free Fermi gas with a 
hard wall at the origin. In this way we obtain 

p{E) = (TtS{E - H)) = 1/5 - siii{2nE/S)/2TTE, (4) 

where S is the level spacing for E ^S> 6. We see that the 
coupling to the superconductor depletes the mean den- 
sity of states and makes it vanish quadratically at E = 
pr| . The states pushed away from E = cause density 
oscillations, which ebb off as 1/E. Note that the result 
applies when both E and S are much smaller than the 
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characteristic energy uncertainty set by the frequency of 
Andreev reflection. 

To understand better the mechanism of depletion, we 
turn to diagrammatic perturbation theory. Let G := 
((E + is — H)^ 1 ) denote the ensemble average of the 
Gorkov Green's function. We expand it in a geometric 
series with respect to Ti as usual. To do the ensemble 
average, we distinguish between two types of contrac- 
tion U acM = v2s ad$bc and U^ c bd = ~v 2 5 aE S bd: corre- 
sponding to the first and second term in the basic law 
(||). Making this distinction is useful for organizing the 
perturbation series, since II D causes pure GUE behav- 
ior whereas II C generates the corrections to the GUE. 
By summing all nested n D self-energy graphs, we get 
Pastur's equation, G — (E + ie — v 2 TrG)~ 1 , which is ex- 
act for E ^> 5 and N — ► oo. The solution, G°, of this 
equation yields Wigner's semicircle law for the density 
of states: -ImTrG /^ = y/2N - (E/v) 2 /irv. Accord- 
ing to (Q), corrections to this result, which is station- 
ary and equal to \/2N /ttv —: 1/6 up to uninteresting 
terms of order 1/N, should appear as we approach zero 
energy. It turns out that these arise from summing a 
geometric series of ladder graphs built solely from II C - 
contractions. The ladder sum, C ac ,bdi satisfies Dyson's 
equation C = C° + C°Tl c C with C° ac bd = 6 ab 6 cd G° aa G° cc . 
Its solution C = C°(l - TlPC )- 1 is singular at E = 0: 

C = C° + C° (Tl c 2iN6/irE) C° + 0(E/N) . 

By evaluating the graph shown in Fig. 3 we get TrG/7r = 
—i/5— 1/2ttE + which are the leading terms in a 1/E 
expansion. (There is a renormalization by a factor of 1/2 
coming from the possibility of connecting the external 
legs in Fig. 3 by a nonsingular II D ladder.) This pertur- 
bative result is to be compared with the exact formula 
TtG/tt = -i/8 - [1 - exp(27ri£ , /6)]/27r£' reconstructed 
from ([|) by causality. We see that diagrammatic per- 
turbation theory properly reproduces the smooth part of 
the 1/E correction. (The oscillatory term is nonanalytic 
in the expansion parameter 1/E and cannot be recovered 
by the perturbative summation of graphs.) 

What is the semiclassical meaning of the mode CI 
Our diagrammatic analysis suggests an interpretation as 
a mode of phase-coherent propagation of a particle-hole 
pair. To gain further insight, recall the periodic orbit 
shown in Fig. 2. In diagrammatic language, this or- 
bit corresponds to connecting the Green's function lines 
in Fig. 3 by just a single Il-contraction, and thus is 
the simplest semiclassical building block of the C-mode. 
More generally, an arbitrary number of Andreev reflec- 
tions may be inserted into the loop. The only condition 
is that the particle-hole character of the states during the 
first and second traversal of the loop be CT -conjugate to 
each other. Because CT-conjugate states carry opposite 
charge and the loop is traversed twice in the same direc- 
tion, the C-mode is ignorant of a weak magnetic field [12| . 



Note, however, that the C-mode is sensitive to energy E 
(or voltage), which breaks the phase relation between 
particles and holes. 




FIG. 3. Diagram contributing to the average Green's func- 
tion in order 1/E. The shaded region symbolizes the C-mode. 

With a solid understanding of the spectral proper- 
ties in hand, we finally open up the tunnel barrier and 
turn to the prime experimental observable, the con- 
ductance g. Our treatment will be based on the lin- 
ear response formula g = 4(e 2 //i)TrS'' i 'pi 1 S'hp @ where 
Sh p is the part of the scattering matrix S that maps 
incoming electrons onto outgoing holes. We suppress 
elastic phase shifts and parameterize the iS-matrix by 
5=1- 2iW{ie + iY- K) _1 W @ at E = 0. Here W 
is a matrix coupling 2M channels (particles and holes) 
in the lead with 2N levels in the quantum dot, W is its 
adjoint, and Y = WW. W is diagonal in particle-hole 
space. For a simple model, we take WW to be a multi- 
ple of the identity in channel space: WW — 7I. Then 
r = 7P where P is a rank-2M projector in level space. 
We assume N 3> M. 

We begin by discussing the limit of an open quantum 
dot (M ^> 1) where all structure in the density of states 
is washed out by the large level width. Let T denote 
the probability for an electron in an incoming channel to 
be transmitted through the tunnel barrier, in this very 
limit. For our simple model, T = 4Xj/(X + 7) 2 with 
A 2 = 2Nv 2 (T4J. Having entered the quantum dot, the 
electron spends a long time there, undergoing many An- 
dreev reflections, and finally returns to the lead either as 
a particle or as a hole, with equal probabilities. Hence 
(TrS^phShp) = MT/2, which we refer to as the "classi- 
cal" value. To confirm this result diagrammatically, one 
has to sum a geometric series of n D ladder graphs, pro- 
ducing the diffuson [^5). On approaching the opposite 
(closed) limit, we expect a reduction relative to the clas- 
sical value. The reason is simply this: the density of 
states in the isolated dot vanishes at E = 0, so that an 
electron trying to enter cannot find any state to go to and 
is immediately rejected into the lead. This effect should 
announce itself as a negative O(M ) correction ("weak 
localization" ) . Such a correction in fact exists and is due 
to the graph shown in Fig. 4 featuring two diffusons and 
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one C-mode. Evaluation of this graph leads to 

(g) = (MT - 1 + T) x 2e 2 /h + <D(1/MT). (5) 

Note that the weak localization correction in the corre- 
sponding N-system vanishes. A result similar to (g) was 
recently found in a related context fll6fl. 




FIG. 4. Diagram contributing to the average conductance 
in order M° . The shading of the vertex center indicates uni- 
tarity-conserving (Hikami box) corrections. 

Although the calculation of the graph of Fig. 4 is some- 
what technical, the result can easily be checked and in- 
terpreted. Consider first the limit of many weakly cou- 
pled channels (T — + 0, MT fixed). In this limit the 
absolute square of the Green's function becomes inde- 
pendent of the level and particle-hole indices, on aver- 
age Q. By conservation of probability, the conduc- 
tance is then fully determined by a one-point function: 
(g) = -(4e 2 /fr)7lm(Tr(i7 _ H)' 1 ) where 7 = jM/N is 
the mean level width. Analytic continuation of our re- 
sult for the average Green's function to E = ij gives 
(g) = (MT - 1 + e~ MT ) x 2e 2 //i, in agreement with (g). 

To understand the other limit of strongly coupled chan- 
nels (T = 1) it is easiest to argue directly at the S- 
matrix level. By its definition as a "sticking probability" , 
T = 1 - |(5 CC )| 2 + 0(1/ MT). Putting T = 1 therefore 
amounts to assuming a fully random S-matrix that van- 
ishes on average: (S) = 0. The symmetries of S are 
the same as those of expiH, so that from (|l|) we de- 
duce S = CS~ lT C~ 1 , which is the defining equation of 
the symplectic group Sp(2M, C). (C is now meant to 
operate in channel space.) Unitarity fixes a subgroup 
Sp(2M). We are thus led to postulate an ^-matrix en- 
semble with probability distribution ((•)) given by the 
Haar (or invariant) measure of Sp(2M). This ensem- 
ble is Dyson's Circular Unitary Ensemble H adapted to 
fit our particlc-holc symmetric situation. By elementary 
group theory, ((S ab )) = and ({S ab S*J) = 5 ac 5 bd /2M, 
so ((g)) = 2Me 2 /h, with no correction of order M°. 



In conclusion, we mention that the singular modes, the 
diffuson and the C-mode, translate into low-energy de- 
grees of freedom of a corresponding nonlinear a model. 
They determine the attractive (metallic, or free) renor- 
malization group fixed point of this field theory. Ulti- 
mately, the existence of such singular modes, which are 
forgetful of microscopic detail and saturate the long-time 
and long-distance physics, is the reason why we predict 
with confidence that the ergodic limit of the chaotic An- 
dreev quantum dot is universal and our random-matrix 
theory applies. 
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